WELL-POSEDNESS OF EINSTEIN'S EQUATION WITH 
RED SHIFT DATA 

CHRISTOPHER J. WINFIELD 

Abstract. We study the solvability of a system of ordinary differential 
equations derived from null geodesies of the LTB metric with data given 
in terms of a so-called redshift parameter. Data is introduced along these 
geodesies by the luminosity distance function. We check our results with 
luminosity distance depending on the cosmological constant and with the 
well-known FRW model. 



Introduction 



Resulting from the Lemaitre-Tolman-Bondi metric 
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are so-called symmetric dust solutions to the Einstein equation given by 

^'■'^ [r) =W + W 

M'{r) 



(0.3) p{t,r) 



R{t,rYR'{t,r) 

(c.f. [17]) for some suitable p (energy density) where superscript / and • denote 
partial derivatives with respect to r and t, respectively. Setting a =^ sgni?, 5 
sgnii', A B^/T+2E and B =^ a^/2E + 2M/R, we study resulting system 
[6, 3] 

(0.4) ^ -^^^ ^ 



dz {l + z)d^^tR{t,r) E' + M'/R- MR'/R^ 

,^ , dt -\R'\ -BR' 5 

(0-5) 77 



dz {l + z)dltR{t,r) E' + M'/R- MR'/R"^' 

taken along null geodesies of (0.1). Here data is given for the function R, pre- 
scribing values R{t{z),r{z)) along curves given by (0.4) and (0.5). As a result, 
corresponding solutions of this system provide maps 

(0.6) iE{r),DLiz),Roir)) ^ (r(z), i(z), M(r(z))) 

as introduced in [6] which we study in some detail in this article. 

As an application of our analysis, we will consider data given in the form 

(0.7) R{t{z), r{z)) = = 

for Dl(z) = {1 + z) Sl''^I{y)dy with I{y) = + (1 - n^)y^ for a real 

parameter < Qa < 1- Here, £)l is generally referred to as "luminosity distance" 
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and, in particular models, Oa = tst is directly proportional to the so-called "cos- 
mological constant" A [3, 5]. For further details of the physical and mathematical 
derivation of the present problem, the author recommends the aforementioned 
articles along with [7, 9, 10, 15, 19, 14] - to name but a few. 

This work is physically motivated by competing cosmological theories in ex- 
plaining certain observations of matter distribution and cosmic inflation. Such 
theories include those of "dark energy" [8] , certain metric perturbations from the 
FRW model [11, 16, 13], radial inhomogeneities of the unperturbed LTB model 
(via E{r), Ro{r) and M{r)), and the cosmological constant (here via Dl) - with 
our work involving the later two. Here, we study the map (0.6) mostly on purely 
mathematical grounds, presenting a framework of analysis and, in a special case, 
estimates on resulting functions M in terms of z. Furthermore, we test our results 
for certain functions Dl, E, and Rq, arising from various FRW- type models, and 
study singularities of M as indications of (in-) compatibility of these models. 



1. Singularities 

From the Chain Rule (c.f. equation (14) [6]) we observe that R' takes the form 



R' = T{R, Ro, R'o, E, E', M, t) + M'g{R, Rq, E, E' , M). 
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With it = a-s/2E + 2M/R and Ro{r) =^ R{r,to) for a fixed to > 0, we restrict 
R, Ro,t>0, M > -ER, E>0 and set 



(1.1) J{R, Ro, M, E, t) V2{t -to) -a £ \j m '^^ 







solutions of which define smooth manifolds depending on constant a = ±1, 
respectively. 

We introduce notation: For a given function / = f(t,r), depending implicitly 
or explicitly on {t,r), we will denote f[z] =^ f{t{z),r{z)) and, with slight abuse 
of notation, set ^ ^=lf ffl. We now set 

' dz az 

where, from the chain rule, with subscript denoting the associated partial deriv- 
ative, 

(1.3) -{dRj)T = E'OeJ + RodR,J 

(1.4) -{dRJ)g = OmJ. 



with C = M/R, ^i, = M/Ro, f) = Rq/R, and 



(1.5) JR = -a/^E + ^ jR,=al^/E + 

Jm — —-;: / TT^ — , ^xo/n Je — —- 
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Substituting (1.2) into equations (0.4) and (0.5), we obtain 



We then substitute 

R' — _ 

R'^" J dz \ R i?2 

so that equation (1.7) becomes 



jr I _ = ^ _ I 



, X .dt , „ ( dr dM ^ 



Equation (1.8) can be verified by equations (0.2) and (0.3). 

Now, equations (1.6), (1.8), and (1.2) along with the Chain Rule result in the 
following system: 

R^-^J dz^\RR^)dz 
,^^dr ,dt ,^^dM 

8BT—+A—+5Bg— =0 

dz dz dz 

■^dr /—— — --TTT-^dt ^dM dR 
T— + asJ{2E + 2M R)— + = — 

dz dz dz dz 

which we may write in matrix form as 

(1.9) «f = y 
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for 



u = 



771/ M 77 





A 



5BF 

T a^/2E + 2M/R 



j_ _ Me 

R 'W 

sgB 
g 



X 



( \ 

r 



dR 
\ dz J 



We check the invertibihty of U as we compute 



M 



detU ={E' - -^J^){Ag - aSgBy/2E + 2M/R) 
+{^ - ^){5gB:F^2E + 2M/R - TA) 



=B{E'g - -){Vl + 2E - a5^j2E + 2M/R) 



Prom these computations we conclude 

Proposition 1.10. Suppose that E,R> Q with dtR, drR, df,.R 0. Then, 

is a smooth function of R, Rq, R'q, E, E', and M except for the following cases: 

Either 

1. ) both 6 = a and R = 2M; or, 

2. ) E'Rg = T. 



We may extend the domain of U to include —1/2 < E < 0, say, but for 
simplicity we impose the above hypothesis throughout the rest of this section. 
We continue with 
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Proposition 1.11. Suppose that for some z* > 0, 5[z*] = (7[z*] and that 
= 0. Then, the matrixU[z\ is singular at z = z*. 

Proof. We have from the Chain Rule and equations (0.4) and (0.5) that 

(1.12) ^ = R'^ + R^ 

az az az 



^ VT+2E ^ \R'\ 



(1 + z)dltR (1 + z)dltR 

_ , V1 + 2E- Sa^y2E + 2M/R 
{l + z)dl,R 

By our hypotheses on the partial derivatives of R we may conclude 

(a5^/2E + 2M/R^ [z*] = (Vl + 2E^ [z*] 

With a = 5 a.i z = z*, we have that crS = = 1 and that ^J2E + 2M/R 
= y/l + 2E so that 2M[z*] = R[z*]. Then from Proposition 1.10 we see that 
det;7[2*]=0. □ 

We note that the type of singularity of item 1) of Proposition 1.10 appears 
analogous to that of the well-known "Schwarzschild" singularity: It is not yet 
clear here if this is merely an artifact of the specific model or if such singularities 
are removable by passing to alternate coordinate systems or metrics (c.f. §31 
[14], §6.4 [18]), taking us beyond the scope of the present article. 

We may interpret item 2) of Proposition 1.10 in terms of the tangent bundles 
TO^ (resp.) of manifolds obtained from (1.1). We may consider the transforma- 
tion 0± : R X (0,+oo) D± given by (/>=^(t,r) =' {R{t,r),Ro{r),E{r),M{r),t) 
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and as a push forward, to interpret corresponding solutions to 

E'ROmJ = E'OeJ + R'odR,J 

as subsets M."^ of T£)^ (resp.) in coordinate form. Let T denote the set (0+ \ 
■kM.'^) U \ ■kM.~) where tt denotes the natural projection tt : TM. M. of 
a manifold A^. 

By calculating U~^Y from (1.9) with '= ^^j^, we arrive at the following 
system of ordinary differential equations: 

dz gE'R-T R ■ {8(Jsj2E + 2M/i? - Vl + ^E){QE'R - 7) 

dt _ 5-R,, 
dz " 5a^2E + 2M/i? - 

dM _ - i?2 • [R^E' - FM)y/l + 2E 

dz~gE'R-T R ■ {8asj2E + 2M/i? - Vl + 2E){QE'R - T) 

We are ready to state 

Proposition 1.14. The matrixU is non-singular for R ^ 2M provided {R, Ro, E, M, 
e T. Indeed, if for some zq > 0, these conditions hold for {R,Ro,E,M,t)[z]\z=zo^ 
then the system of equations (1-13) has a unique C°° solution X[z\ in some open 
interval containing zq. 



Proof. It is clear that the elements of U are continuously differentiable where 
detf/ is non-zero. The result follows by applying standard theory of ordinary 
differential equations [4]. □ 
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To further investigate the solvability of the system (1.13), we compute 



(1.15) GE'R-T- 



E'R ■ [Je/ R — Jm) + RqJrq 
Jr 

r^r—jRE^ s^l\s-l)ds % 



Since jeU^ ^ 27W 



1 {Ei^ + C)3/2 - ^V27^ Ji 



Lacking any other simplifying assumptions, we thus obtain strong criteria for 
local solvability: 

Proposition 1.16. System (1-13) is locally solvable at any point of where 
S ^ a or where 2M ^ R if either of the following holds: 



1) sgnE' ^ squRq 

2) I 



Indeed, given ro,to,Mo > and smooth E, R, R'q > 0, the system (1-13) has on 
an open interval I B zq a unique solution satisfying 



X{zo) 



( \ 

to 

y Mo y 
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2. Decoupled equations: A Case of Constant E 

We consider the case of constant ii^ > in which we can rescale M and R to 
assume the case E = 1, retaining 



Here, equations (1.13) reduce to 



(2-2) :r = ^^ + 



dr -GAR , i?^ ■ (1 - f g)x/3 



dz T (V3 - (j5y/2 + 2M/R)T 

dt _ -()■ 

dz~ s/^- ab^2 + 2M/i? 

dM , Rz^^ 



= RA 



dz ^ - a5^/2 + 2M/R 

withA=^^^^#^. 

l+z 

For the remainder of the section we assume that E,a,S = 1 and denote by 7i 
the corresponding subset of T. Then, R{t,r) < Ro{r) Mt <to. And, for f) and ^ 
as above, we obtain 



OmJ = 
dRaJ = 



1 

2. /I 



1 



1 



1 + e 
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with C > and f) > 1, so that the following hold: 



(2.3) o<Mr,z,o = -j = ^^l iVTlVTl'" 

^ V^f)(l + ^«) / 1 Vi) . 



Our change of variables leads to 

dM 



dz ../^-^^«/^ 



with A = ^^^^ whereby the system (1.13) now reduces further to 



(2.4) 



dz 1 + z V3 - + 

dt _ -R^ 
~ V3 - V2 + 2^ 



dz 1 + z i? \^y3- v'2 + 2^ J ' 

Here, we note that the equation for ^ decouples from the others, allowing for 
^ to be solved for explicitly in z. Then, with the solution to ^(z) in hand, both 
Ii and I2 depend only on z and r whereby the remaining equations are then 
decoupled. 

We give estimates for the system (2.4) assuming uniform bounds on R, M/R, 
Rz, Ro, and R'q. We suppose the following bounds hold for < zq < z < Zi and 
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< r,M,t on some compact sets (to be determined): ^ < ^* with |2^ — 1| > e 
> 0; Pmin < R< Pmax; \Rz\ < A ; 1< f) < 1^*; and, > r > 0. Here, applying 



(2.3) 



(2.5) 



dz' - ^3 - V2 + 2^ 



A-(V3 + V2 + 2e) def 



yr72 + (1 + yF/2)mti def 

1^1 < \/3(i+e)(V2+$Hir/p™„) ='wt3 



Let 9}t maxj{9}tj}j^i and suppose tq, to, and Mo/-R[2;o] =^ Co 7^ 1/2 satisfy 
the restrictions on (r, C) for some < ^o < C* as above with 



(2.6) 



^0 = X(2o) 



to 



For an interval / of the form Q < Zq < z < zx, the following now results from 
standard theory of differential equations [4]: 



Proposition 2.7. For zq > 0, the system (2.2) is solvable on an interval of the 
form, I = {z\zo < z < zi} provided that the conditions (2.6) and (2.5) hold for 
X in subset of 71 given by \{X — Xo)j\ < b : j = 1,2,3 for some constant b 
< 1/9JI. Here, a unique solution may be computed by the method of successive 
approximations. 
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Proof. We may apply Theorem 3.1, Chapt. 1 [4]: The conditions assure Lipshitz 
continuity of the right-hand sides of (2.4) and that both z — and \X — Xo|/SOT 
are bounded above by jzi — Zo|, so that the result follows. □ 

Recalling that we set = 1, we will suppose for the rest of the section that 
R[z\, R'q[z\ are smooth and positive for z> Q. For some of our analysis below we 
will suppose also that 

(2.8) R[z] > Cz\R^\. 

holds on some real interval. We now present our estimates on M[z] depending 
on R[z\ and initial conditions given by '= ^('"(-^o)) -^o)- 

Theorem 1. Suppose that (2.8) holds on some interval I = [zq,zx) C M"*". 
Then the following statements hold for some constants 0<ci<l/2<C2, each 
depending on the choice of C: 

1) 7/0 < ^0 < 1/2 and < 0, then M[z] < ciR[z\ holds on I. 

2) 7/^0 > 1/2 and Rz > 0, then M[z] < C2R[z] on I. 

Proof. Let us choose C < 1/2 and set A| =^ \/3 — \/2 + 2^. In case 1) we use 
the estimate l/Aj > ^3/(1 - 2^) for < ^ < 1/2 so that from (2.4) 

Here, f < for 1/2 > ^ > ^Jf C/{2C + 1). Let ci = max{$o,Ci*}. 
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In case 2) we note that 1/A| < -\/3/(2^ - 1) for ^ > 1/2. We find 
and f < for 1/2 < ^ < ^2 C/(2C - 1). Let C2 = maxl^o, □ 



Theorem 2. Suppose Rz < on I = [-20,00) with zq > and ^0 > 1/2. T/ien, 
for p =^ -^3/2, i/iere are positive constants a, C3 and C4 so i/iai i/ie following 
holds on I : 



* * < < C4 ' ^ ' 



+ zoJJ- "-*\{R[z])(-P-y^') 

Proof. We first note that since Rz/A(^ > on /, wc find from (2.4) that ^ > 0. 
Let us set =^ -y/Sy/T+J/A^, noting that p =^ y/3/2 < < for ^ > 1/2 
is decreasing as function of ^ and, in turn, also as a function of z. Recalling that 
Rz < 0, we find 



d£ Rz Rz /^ VTT^ 

dz R[z\ dz R[z\ 1 + z 



Now, by the monotonicity of R[z], 



.0 - + « 



> R-P[z] (^^oR'[zo] + RP[z]V6^/T+To j^' 



After multiplying through by R, it is clear that we may choose C3 < min{^o(-R[-zo])'', V^Vl + ^o}-| 
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Now, let us set po K^^ and qo =^ (^o + 1) /Co- Then, for obvious substitution 



defining ^[s] 



dz R[z\ 1 + z 



< R-P^iz] (^^oR'°[zo] + RP"[zo]V6V^n[z\ J^^ Y 



ds ^ 



CM < qoR^f"'[zo]R-^f"'[z] (l + VeJ^^ - ^ ^ 

noting that C^/(l + C) — C/?o- Choosing C4 > qo6R'^P°[zo\, the result follows by 
multiplying through by R. □ 

We see that Theorems 1 and 2 can apply for R[z] = Rq^ [z] (modulo a rescal- 
ing factor) as above for certain values of Qa- We denote by 7^^ the subset of 
(0, 00) for which ±Rz > and we replace C by in the case that (2.8) holds, 
respectively. 

Remark 2.9. For -RaAf-^] ^ (0-7) we find that when 0,^ = 1 there is to every 
interval of the form (0, Z2), an associated C+ depending on Z2 > 0. Moreover, for 
every < f^A < 1 there is a za > where for every positive z"^ with z^ ^ za 



16 



CHRISTOPHER J. WINFIELD 



there is a associated to (0,2;+) and {z ,00), respectively. [The singularities 
z\ will be discussed in further detail in Section 3.] 

Proof. For Oa = 1 we find = (0, 00) with zRJR[z] = l/{z + l). For < f^A < 
1, it is not difficult to show that z/R[z] is bounded from below on (0, 00) by a 
positive constant, depending f^A- Therefore, the sign of zRz/R[z] depends on that 
of Rz. For Oa < 1 we find that the sign of Rz is same as that of {z + 1)1 {z + 1)— 
J^^^X{y)dy which is a monotonically decreasing function of z with a unique 
positive root za > 0, depending on Ha- So, Iq^ = (0, za) and /^^ = {zatOo). 
Hence, for 2;+ as above, there are positive constants C+ so that zRz{z) / R[z\ > 
±C+ on intervals (0,2;+) and (2;", 00), respectively. □ 

In a certain case of interest, we find that for certain initial conditions the 
growth of M[z] roughly follows that of a power function for large z. 

Corollary 2.10. In the case of Theorem 2 we have for R = Ra^ with < 
VLa < 1 that, given Mq > 2i?[2o] > and zq > za, for any a > there are 
positive constants ki and k2 so that for p = \^S/2, 




kizP-^ < M[z] < k2z'^P-^+' 



on I = [20) 00). 
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Proof. It is not difficult to show that to any such Oa there are positive constants 
Ci and C2 so that 

Ci/z<IkiM<C2lz 
holds on I. The result immediately follows by Theorem 2. □ 
We may also conclude 

Corollary 2.11. If either case 1) or 2) of Theorem 1 holds on I = [zq,z{), then 
r{z) and t{z) are both solvable on I. Moreover, r{z) is strictly increasing and t{z) 
is strictly decreasing on I. 

Proof We find that and are smooth functions of z since 1/2 is smooth. 
By inspection, we find that ^ < on J so that, by our assumption on a, we see 
for f) as in (2.3) that f) > 1, increasing with z. Then, I2 + ^Ii > Q for z & I and, 
hence, from (2.4) we see that ^ > for ^ e /. □ 

We note finally that these results are consistent with physical interpretation 
where t is interpreted as "look-back" time from an observer at r = with a 
(locally) expanding universe (c.f. [5, 12]). 

3. Study of Singularities, part A: Critical points depending on Oa 

We now consider how singularities may depend on the parameter Qa for R{z\ 
= Rqj^ [z] . As in Proposition 1.11, a singularity arises at ^ = za where 

rpi (l+ZA)-X(l+ZA)-/r^"X(t/)dt/ 
(3.1) [R,\^^^^ = ^^^^-p = 0. 
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Proposition 3.2. The values z\ satisfy z\ > 1.25, increasing as a continuous 
function o/Qa 'in the domain < Oa < 1- Moreover, there are positive constants 
Ci , C2 and C3 so that 



1 V4 1 
< ZA + 1 < C3 



V Oa- Hence, za — +00 as Cl\ 1. 

Proof. It is not difBcult to show from (3.1) that -zaU^^o = 1-25 and that z\ > 
VOa- Now, let us set g =^ 1 + ^a and note that (3.1) gives q2{q) = J^2(y)dy. 
Implicit differentiation now gives 



d^A dq Ji 90a d^A 



Applying ql^{q) =I^{q) I{y)dy on the second term, right-hand side, we com- 
pute 

where K{y,q) = X^{q){q^ - 1) -l^{y){y^ - 1) and Q{q) ='' 3q^{l - Oa)- Since 
I^{y){y^ — 1) is strictly increasing as a function of y > 1, we find K{y,q) > 
for 1 < y < g. Thus, ^ > VOa and hence q > 2.25 VOa- For fci 
4 j"^ "^^ 2(ii)K{y, 2.25)dy we estimate 

so that 
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and our choices of ci and C2 are clear since 



> fci In ( ) + 2.25* 



Next we note that 

dy 



so that y/q < ^-^_^q^ with 1/^2 =' J^'^^ ^ . We choose C3 = A:| and we are 
done. □ 



With a broad range of values ^Qa) bound by the estimates of Proposition 
3.2, one may expect difficulties in applying the present work to cosmological 
models - with singularities zqj^ well within observed redshift values [3, 5, 6]. 
However, some such singularities may conceivably be of type 0/0 if both 
and + 2E[z] — ^/2M[z]/R[z] + 2E[z] were to have zeros of identical order, 
rendering the singularities, in some sense, removable. We demonstrate such a 
case in the next section. 

4. Study of Singularities, part B: FRW Model 

Using solutions from the well-known Freedman-Robertson- Walker model, we 
analyze our map {E,Dl,Ro) — >■ {r,t,M) and study singularities of the system 
(1.13) and their dependence on SI a- We restrict the map as follows: We fix the 
function E{r) and restrict Dl{z) and i?o(f') to certain one-parameter classes in 
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the pre-image space; and, we fix the function M(r) in the image space. Here, we 
consider data given by R[z] = Rnf,[z] as in (0.7) and we set 



(4.1) E=-,M=-,Ro = cM/E = cr 



for parameter c > 0. Following [1], we have R{r,t) = r-a{t) where for some (real) 
parameter r] with kc '= Vc + c^, 



cosh 71 1 dcf 

(4.2) a{t) = h (ccosh?7 + kc sinhrj) = ddl) 

t=— j= — ^ + V^(c sinh 7/ + fcc cosh ?7) =^ (Sg(jj) 



Here, rj is known as "conformal time" which in our case depends on a and thy rj = 
ri(t) = ^a(T) ' &&c^ invertible for rj on an open 

interval containing 0. In particular, is invertible for r} > — arctanh(2fcc/ (l+2c)) 
and C5c is invertible where a > 0, so that a{t) =^c° ®J^(i) indeed holds for t in 
a neighborhood containing kc- Moreover, using (0.5) and setting c = o(fo) with 
to '= t{zo) = V^kc for some zq > 0, 
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Given R[z], we find, indirectly, the resulting solutions of (1.13): 



(4.3) 



t{z) = <3,{d-\a[z])) 



^{z) = R[z]la[z] 



(1 + Zo)C 



M[z\ 



1 / jl 'I{y)dy' 



2 \^ (1 + zo)c J 

As for the relevance of this case to physical models, we note that the associated 
energy density p[z] is a smooth function on (0, oo). 
We are ready to state 

Theorem 3. For any given < Oa < 1 and 2;o > there is a smooth function 
Ro{r) so that E{r), M{r) as in (4-1) and R[z] = Rqj^[z], the system (1-13) with 
initial conditions 



(4.4) 



^(^o) 



R[zo]/c 
V2kc 
\{R[zo\/cf/2j 



has a smooth solution X on an open interval I 3 zq. 

Proof. For those r] where the solutions (4.2) hold we also have R' = a{t) > and 
it' = <^ = > 0. Since the initial conditions hold for r/ = 5-^(0(^0)) = 0, 

(4.3) also holds for r} in some interval containing 0. From continuity arguments 
we see there is also some open interval I 3 zq on which such solutions X{z) in 
turn hold. □ 
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We note that the above method provides no solutions for r{z) and M[z] in the 
case -Ro = unless more data is prescribed, such as asymptotic conditions for the 
ratio R/c in terms z and Zq (c.f. Example A, p. 5 [6]). Moreover, we note that 
singularities may occur in the form R = 2M and/or d = away from zq so that 
we may not arbitrarily extend the domain I of the solution via Proposition 1.14. 

We may apply Proposition 1.14 in regards to uniqueness of solution: To rule 
out one type of singularity, we compute E'QR —T via (1.15). First, we set 
^ = r/{2a{t)) and = I'/i'^c) and compute 

E' ■ [Je - RJm) rE^RE' T^"^*^ «'^'(^ - ^)ds 



Jr ^ ' ^ 2 7i {Es + 0'^^ 



1 r v^(-«w),.,m±i<o 



SJaW (r + 1)3/2 y a{t) 
for c, a{t) > 0. We now compute. 



RqJro / Ro I / R 



Jr V £^^0 + M' \ ER + M 



\a{t) 



c+iy a{t) 

which is strictly negative. Therefore, E'QR —T < and we have ruled out case 
2) of Proposition 1.10. Knowing also that ^'[2;]U=2o 7^ in this case we state 

Theorem 4. The solutions of Theorem 3 are unique for zq ^ z\. 



The solutions (1-13) stand in glaring contrast to the result of Proposition 1.10: 
Indeed, we note that the right-hand sides of equations (0.4) and (0.5) under the 
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conditions of Theorem 3 have no positive singularities zk as £" +M' /R —MR' / R^ 
= r + r/a > 0; yet, we find that the determinant of U in (1.9) vanishes at 
z = za- Since Theorem 3 apphes in the case Zq = za, one may suspect that these 
singularities are, in some sense, removable - so we shall see in remainder of this 
section. 

We give specific cases, depending on Rq, in which the solutions X can be 
smoothly extended across singularities z = za- For such solutions to be valid, it 
suffices that d[z] > is smooth, that (3.1) holds, and that as in (1.12) R[z\] = 
2M[za] (or perhaps as smooth extensions defined at za)- Then, 




((l + ^o)c)^ 
(1 + zaT 



and, hence. 



J(l + za) = R[za] = 



((l + 0o)c)3/^ 

(1 + ^a)3/2 • 



From this we obtain the corresponding value of c by which we define 



(4.5) 



dcf 
CA = 



1 + .--A 



(1 + ^o)(^^A + (1 - ^^a)(1 + ^a)3)V3 • 



We are ready to state 



Theorem 5. Under the hypothesis of Theorem 3 for every > and < CI a 
< 1 there is a smooth Ro{r) for which the resulting solution X{z) with initial 
conditions (4-4) can be uniquely extended to be of class C"((0,oo)). 
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Proof. We take -Ro('') = c^r for as in (4.5). Using (4.2) and following the 
Chain Rule formula 



(4.6) a{t{z)) 



da\z\ 
dz 



da\z\ 



da[z\ 
dz 



it suffices to show that r][z] is smooth and that ^ is strictly positive on (0,00). 
To do this, we set 

with dt[z] =^ ^^^dz, and we proceed to analyze the integral. We may write 



(4.7) 



r{z){r{z)'^ — a{z)) 



for some real- valued function H > 0, analytic for 2 > 0. Here — o is an 
increasing function which vanishes at z\ and is of the same sign as that of —R^ 
y z>0. 

Now, we check the behavior of ^ near the singularity, applying analyticity 
arguments as follows: Using (4.3) and (3.1) we compute 



d^R[z] 



dz^ 

dM[z] 



3(1 -Oa) 



(1 + ^a)2:^(1 + ^a) <0 



dz 



and, in turn, we find that 



3J(1 + ^a) def 
2(1 + za) 



fSlA > 0; 



2M[z] - R[z] = 29JtA • {z - za) + ^i{z){z - za^ 
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R\z] 

for some analytic function Here, ' < on a neighborhood of z\ where 
Rz has a zero of order exactly 1. We therefore find that the following limit exists 
as we compute: 



dR 

lim ^ 



Z->ZA 



2M[z] - R[z] 



d'^Rjz] 
dz-' 



o dM[z] 
dz 

(1-Qa)(1 + ^a)X(1 + ^a) 



> 



We may conclude therefore that 2m[z]^r[z] extends to an analytic, positive- valued 
function on (0, oo) and, hence, C"((0,oo)) 9 > 0. Therefore, a[z] is well- 

defined and is non-zero; and, moreover, r][z] is of class C"((0, oo)). The unique- 
ness follows since Theorem 3 applies to any open interval not containing z\ □ 

Remark 4.8. Following the calculations in the proof of Theorem 5, we note that 
any other choice of positive c ^ ca leads to a singularity of order one at z = za 
for R'[z] = d[z] as evident in (4.6) and (4.7). This gives singularities in equations 
(0.4), (0.5), and (1.2), and renders the resulting system (1.13) invalid at za- 

Remark 4.9. Our FRW model is consistent with the construction of R[z] as in [5] 
where ^ = 1 + z with no prescribed value of c. Moreover, our choice of c = ca 
is optimal in assuring the largest possible domain of C"^ -solvability. 

Discussion 

We make several concluding comments and a conjecture: First, we note that in 
the case of Theorem 5 the various right-hand sides of the system (1.13) can each 
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be written in the form ji^z]-2Miz] smooth functions A and B. Thus, 

the arguments for the smooth extension of ^ beyond the critical points z\ of 
R[z] also apply to ^ (also, of course with circular reasoning, to However, 
in the general mapping scheme (0.6), we have no way to predict the order of 
the zeros of ^ nor any a priori justification to expect these singularities to be 
removable - not even as we fix our choice of R[z] = i?OA[^]- 

Second, one may interpret the removability or non-existence of such singular- 
ities as indication of compatibility of the corresponding models as one imposes 
R[z] on a model that prescribes E and R^. (Here the LTB model would be said 
to 'mimic' the given cosmological constant model, c.f. [2].) Applying such crite- 
ria to Remark 4.8 one does not expect every LTB model to be compatible with 
such a cosmological-constant model (at least not for z near za)- However, from 
Theorem 5 we do find, as a check of our analysis, that the cosmological-constant 
models for < IIa < 1 are each compatible with at least one LTB/FRW model: 
Our choice of Rq identifies an optimal FRW model, in the sense of Remark 4.9. 

Finally, one conjectures that these removable, 0/0-type singularities may yet 
lead to instability of numerical solutions of the system (1.13) (but here at certain 
finite z (!) c.f. §IV [6]). Such investigations are beyond the scope of the present 
work. 

References 

[1] H. Alnes, M. Amarzguioui, O. Gron, Phys. Rev., D73, 083519, 2006. 



WELL-POSEDNESS OF EINSTEIN'S EQUATION 



27 



A. Aguirre, Z. Haiman, Cosmological Constant or Intergalactic Dust? Constraints 

from the Cosmic Far-Infrared Background, ApJ, 532:28-36, 2000. 

M.N. Celerier, Do we really see a cosmological constant in supernova data?, A&A, 2, 

2008. 

E. Coddington and N. Levinson, Theory of Ordinary Differential Equations, McGraw- 
Hill, New York, 1955. 

S. Carroll, W. Press, E. Turner, The cosmological constant, ARAA, 30, 1992. 

D. Chung, A. Romano, Mapping luminoscity-redshift relationship to LTB cosmology, 
Phys. Rev., D74, 10103507, 2006. 

K. Enqvist, Lemaitre-Tolman-Bondi model. Gen. Rel. Grav. 40:451-466, 2008. 

J. Preiman, Lectures on Dark Energy and Cosmic Acceleration, AIP Conf. Proc, vol. 

1057, pp. 87-124, 2008. 

J. Islam, An Introduction to Mathematical Cosmology, Cambrige Univeristy Press, 
2002. 

J. Kristian, R. Sachs, Observations in cosmology, ApJ, 143, 1966. 

E. W. Kolb, S. Matarrese, A. Riotto, On cosmic acceleration without dark energy, 
NewJ.Phys. 8:322, 2006. 

E. Kolb, M. Turner, The Early Universe, Addison Westly, 1990. 

J.W. Mofit, Late-time Inhomogeneity and Acceleration Without Dark Energy, JCAP, 

0605, 001 (2006) 

Gravitation, C. Misner, K. Thorne, J. Wheeler, W.H. Freeman and Company, 1973. 
H. Pajrtovi, B. Mashhoon, Toward verification of large-scale homogeneity in cosmology, 
ApJ 276, 1984. 

S. Rasanen, Backreaction in the Lemaitre-Tolman-Bondi model, JCAP, 0411:010, 
2004 



CHRISTOPHER J. WINFIELD 



] H. Stepani, D. Kramer, M. MacCalUum, C. Hoenselaers, E. Herlt, Exact Solutions of 

Einstein's Field Equations, Cambrige University Press, 2003. 
] R. Wald, General Relativiey, University of Chicago Press, 1984. 
] S. Weinberg, Gravitation and Cosmology, Wiley and Sons, Inc., 1972. 

E-mail address: winfieldOmadscitech.org 



